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Perturbation theory and harmonic gauge 
propagation in general relativity, a particular 

example 

Jesiis Martin* Alfred Molina^and Eduardo Ruiz^ 


Abstract 

We study how the changes of coordinates between the class of harmonic 
coordinates affect the analitycal solutions of Einstein’s equations and we ap¬ 
ply it to an analytical approach for stationary and axisymmetric solutions 
of Einstein equation used by the authors m. m to solve the problem of a 
self-gravitating rigidly rotating perfect fluid compact source. 


1 Introduction 

In several previous articles BE] to solve Einstein’s equations we used a new method 
to obtain approximate stationary and axisymmetric solutions to describe the grav¬ 
itational field (interior and exterior) of a compact stellar object in rigid rotation. 
This method is a combination of a postminkowskian approach and another one that 
could be qualified as post spherical. At every postminkowskian order in a dimen¬ 
sionless parameter, A that we set as usual to be proportional to the mass of the 
object, a series expansion is performed in powers of a parameter, C, also dimen¬ 
sionless, which takes into account of the deformation of the object relative to the 
spherical symmetry. This parameter classically represent the ratio in the equato¬ 
rial plane between the centrifugal force and the gravitational one, or in terms of 
the Maclaurin ellipsoids some function of the eccentricity. In this sense this 

second subordinate series expansion is similar to the slow rotation approach used 
by Hartle |S], but with the advantage that we do not start from an spherical exact 
solution in A and then we can go easily to higher hi terms and we do not need any 
numerical computation to obtain algebraic results in the slow rotation parameter. 

In this article we are going to study only some aspects related to the post¬ 
minkowskian approach without any care with respect to the post spherical approach. 
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To explain our current goal let’s recall first that the approximate general solution 
of the Einstein equations to the postminkowskian linear order in harmonic coordi¬ 
nates contains four groups of constants. As regards the exterior solution two of 
these groups of constants are the static and dynamic multipole moments of Geroch- 
Hansen, Thorne et al P, [71 [H] (and therefore they are invariant), while the other 
two groups are removable by a coordinate change, so that they can be described by 
gauge constants. On the interior a similar situation occurs, there are two sets of 
non removable constants, but they do not have a specific meaning as for the exterior 
solution, and two other groups gauge constants. A priori it seems natural to choose 
a coordinate system in which all gauge constants become equal to zero, this is the 
choice of ThornejS] and that is known as canonical gauge. However we are interested 
in making a match on the surface of the stellar object by using the Lichnerowicz’s 
prescription [10], i.e. in this matching on the surface the metric is continuous and 
has continuous derivatives. On this basis the gauge constants (both external and 
internal) are needed and also are uniquely fixed when the matching is performed [2]. 

As a consequence of this we need to maintain the gauge constants at the first 
order and drag them to the second order (and so on) in an iterative process that 
suppose the postminkowskian approach. It is easy to see that this result in quite 
long and cumbersome calculations, even with the use of the algebraic software, which 
must be implemented with elaborate routines when we go over the second order. The 
following questions can be naturally raised: 

• What is the aspect of the solution at the second order (for instance) that comes 
from the first order gauge constants? 

• You may obtain this part of the second order solution just analyzing inhnites- 
imal coordinate changes up to the second order? 

The main objective of this paper is to answer these questions, showing that 
the global solution to the second order contains the first order gauge constants 
transferred to the second order solution that matches exactly with we expect from 
the analysis of inhnitesimal coordinate changes. One might object that the result 
is trivial because of the covariance of the theory; however this argument present 
some difficulties which are specified in the text, and in any case our contribution 
is a explicit verihcation of how gauge constants are transferred, which can simplify 
the calculations significantly. 

The paper is organized as follows. In Section 2 the expressions for the harmonic- 
ity conditions and the Ricci tensor are written in an appropriate way to develop 
the postminkowskian approximation. In Section 3, these expressions are used to 
write Einstein’s equations to different orders in the postminkowskian development. 
In Section 4 the infinitesimal coordinate changes are studied up to the second or¬ 
der including the peculiarities of the harmonic changes , i.e., those who go from 
harmonic coordinates to harmonic coordinates, specially the results are written for 
our problem i.e. for an spacetime stationary and axisymmetric with a Papapetrou 
structure and asymptotically flat. In Section 5 the properties of the stationary and 
axisymmetric metrics are specified. In Section 6 we address the problem of exter¬ 
nal vacuum solution, solving Einstein’s equations to first and second order of the 
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postminkowskian approximation in harmonic coordinates. Besides we prove the fun¬ 
damental thesis of this work, i.e., the part of the second order solution transferred 
by the hrst order gauge constants can be obtained by a simple analysis of the in- 
hnitesimal changes. It is also argued that the result must be true when we are going 
to the third order. In Section 7 the same process is repeated for the case of an 
interior perfect fluid solution in a rigid rotation and with a linear equation of state. 
Finally in Section 8 the conclusions of the work are presented. 

As usual the space-time metric is written as follow 

ds"^ = ga/3{x^)dx°‘dx'^ 

The Greek and Latin index take the values 

a,/3,A,... = 0,1,2,3 ; = 1,2,3 

We use the Einstein summation convention and the signature of the space-time is 
(“) +) +) +)• 

2 Harmonicity conditions and Ricci tensor 

2.1 Harmonicity conditions 


It is well known that the harmonic coordinates condition is 



r" := = -^9^ (^J<r) = 0 

(1) 

but 

Kp = log 

(2) 

and then we have 

r“ = 

( 3 ) 

We can dehne also 

Ta := = -gaxdf.g'"^ - 

( 4 ) 

If we dehne the metric deviation from Minkowski metric as 




( 5 ) 

and the inverse metric 



^a/3 ^ ^ ^ + k’^^kp^ = O) 

(6) 

the equation (jl]) can 

be written in terms of the metric deviation as follows 



To = a,\p + k^^T = G -|- Pq, = 0 

(7) 

where 



Pa-- = 

- dah 

(8) 


by splitting the linear terms in the deviation from the nonlinear terms. 
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2.2 Ricci tensor 


Now we are going to write the Ricci tensor in terms of the left hand side of har- 
monicity condition and the metric deviation. First of all from the dehnition of the 
Ricci tensor we have 


Rcr> = ajry - (9) 

but using ([3]) the hrst term of the right hand can be written 

= ( 10 ) 

and using this result in ([2]) 

( 11 ) 

In this expression only the hrst and the second term in the right hand side contain 
linear terms in the deviation. We are going to split them from the non linear ones; 
since 

g^^dxT^^o.13 - + 5 a /3 Va - dx^^Kfi) 

+ k^^dxTf,,ai3 - ^^g^^dp^hxf, - (12) 

then the Ricci tensor can be expressed as 


Rap — Lap — + Nap 


where we have dehned 
Lap ■ 


' i^—^Lap ~\~ ^da^p ^dpla 


( 13 ) 


( 14 ) 


N^0 ■= - la^(fc^'‘a„AAj - (15) 

and where all the linear terms in the deviation are in Lap- If we use the harmonicity 
condition Fq = 0 we can substitute la = —Pa and we hnally get 

Rap = --jjLhap — -daPp “ -R^pPa + ^ap (16) 

where only the hrst term in the right hand side contains linear terms in the deviation. 


3 Einstein equations and series development 

The Einstein equations read 

Rap =Tap -^Tgap ■.= TaP (SttG = C = 1) (17) 
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Let us assume that the deviation of the metric and Tais can be developed in a 
series of a parameter A, i.e. 

hap = ^h^al + -= X] ^""h^a} 

n=l 

< 

n- _ X'T'd) _i_ \2q-{‘i) I , _ \ ^ \nq-{n) 

tap - ^tap + ^ 'ap + ^ > aP -- ' aP 

^ n=l 

Then the linear term as in the harmonic coordinates condition ([7]) can be written 

L = Y, := , a<"' := (19) 

n=l n=l 

For the nonlinear term we have 



Pa = >?P^a^ + + . . . = ^ A"PW •= 

n=2 n=2 r-\-s=n 

{r,s = 1,2,...) where in only the terms are needed. The terms 

of the inverse metric are obtained as follows; from the dehnition 


k<^f^ = X^k^^P^ and the relation = 5^ 

n=l 


we obtain 


n=l 




r-\-s=n 


= 0 


If we know this gives us k^^P'" and if we know 
{hW,...,hW} and 


the term with A" give us the term k^P". Then to obtain the term we need the 
deviation of the metric to order n — 1. 

The same happen with the Ricci tensor ffTBD . The remaining terms to analyze 
are in Nap, in which the two hrst terms have the same structure we have already 
analyzed, and the third one contains products of the type 


{v^^ + kP)Tp,ap{v''‘' + knr.,px 


The terms of order n are obtained with the deviation of the metric to order less 
than order n, and like Pa in Nap the lower term is of order A^, i.e. 

+ AWf; + ■ ■ ■ = \‘Np^ (20) 

n=2 


Then to first order the Einstein’s equation and the harmonicity conditions are 


□ — —2 
'''a /3 “ 'aP 

(9'’hpa — - d°‘hP> = 0 


( 21 ) 


5 






and to order n > 2 


□ /ii”j = -2r7 + 2 


(n) 


N. 


(n) 




a pin) 
O(a-r^) 


d^h 


(n) 

pa 


Qa^in) ^ _pin) 


( 22 ) 


2 

'(2) a/-(2) 


where for instance for n = 2 and are build with the linear deviation 

i-e- 


p(2) _ L(l)AAip(l) 

1 1 (23) 

ivi5 = ^ 

and is build from The same happen with the higher order equations, the 
second member of the equations are build with lower orders that are already known. 

4 Infinitesimal coordinates change to second or¬ 
der 

Now we are going to study how the infinitesimal changes of coordinates affect the 
deviation of the metric of the previous section and how this gauge propagates to 
higher order. Let us assume that we take the deviation to second order in the 
parameter, i.e. 

ga^{x) = r]afi + Xha 0 {x) + X^qafi{x) (24) 

and the inverse metric 

g^^x) = + Ar^(x) + AV^(x) (25) 

Let us point out that we have changed the notation, the terms hap and k°‘^ represent 
from now on the first order in the parameter A for the deviation and not all the 
deviation as in equations ([5]) and ([HD. i.e. 

f ^a/3 ^a/3 ) ^a/3 ^ Qa/S 

I ^(2)a/3 pa0 


Then from ([HD 


= 0 ; = 0 


If we perform now the inhnitesimal coordinates change to second order 


x“ = x“ + A^“(i) + A^C“(i) 


( 26 ) 


( 27 ) 
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the metric in the new coordinates can be written as 
go,fi{x) = VafS + + >?qafi{x) 

= Vafi + >^Kp{x) + C^{x)dphafi{x) + g«/j(x) 

but also 

n r,i nr] — \ i 

^ dx^ 


So/jW = 


(28) 

(29) 

(30) 


and then 

^0/3 + ^hctfi^x) + A^5q,^(x) = + \da^^ + A^daC^)(^0 + A9/3^^ + 

X |?7am + ^h^(x) + A^ ^^(x)dphxp(x) + gxp(x) | (31) 

Now we can identify the two hrst orders in A in equations fl2^ and fl3T]) 

ha^(x) = Kfi{x) - daifi{x) - 9/3^0(x) (32) 


qap{x) = qaf}{x) - daCfi{x) - dfiCaix) - dai^{x)df}ip{x) 
- ^^{x)dp hafs{x) - da^fj{x) - dfsCaix) 

haf^{x) - da^p{x) - dp^a{x) dfj^^{x) 

hpfi{x) - dpip{x) - dp^fiix)] dc^f^ix) 


(33) 


which can be written as 


qaf}{x) = qap{x) - daC/3(x) - d^Ca(x) - £(^)/i„/3(x) - da^^(x)d^^p(x) (34) 


where the two last terms are the hrst order gauge propagation to second order. 

To use it later on it is important to specify the equations that should verify 
the vectors ^“(x) and C^lx) to accomplish that the change (31) transforms har¬ 
monic coordinates to harmonic coordinates (“harmonics changes”). As is known it 
is necessary that 

^"^aA^x“(x) = 0 (35) 

which results in the following equations 


v"^dxpC(x) = 0 
rf^'^d;,pC(x) = h^^(x)d;,pC(x) 


(36) 


having taken into account (29), (30) and (31), and where = ri^°‘ri^^^hai 3 ■ 
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5 Stationary and axisymmetric metrics 

We require space-time to be a stationary and axisymmetric Riemannian asymptoti¬ 
cally flat manifold admitting a global system of spherical-like coordinates {t, r, 6, 9 ?} 
which verifies the following properties: 

A) Coordinates are adapted to the space-time symmetry, that is to say, ^ = dt and 
C = dip are respectively the timelike and spacelike Killing vectors, so that the metric 
components do not depend on coordinates t or (f. 

B) Coordinates {r, 6 *} parametrize two dimensional surfaces orthogonal to the orbits 
of the symmetry group, that is the metric tensor has Papapetrou structure, 

-h 7 ^ 0 (n;'’< 8 ) 0 ;® , (37) 

where cj* = dt, = dr, = rdd, = rsm.6 dip is the Euclidean orthonormal co 
basis associated to these coordinates. 

C) Coordinates {t,x = rsin^cos:^, y = r sin 6 * sin (^, z = cos0} associated with the 
spherical-like coordinates are Cartesian coordinates at spacelike infinity, that is the 
metric in these coordinates tends to the Minkowsky metric in standard Cartesian 
coordinates for large values of the coordinate r. 

All these properties are compatible with the harmonic coordinates we use in this 
paper. Time coordinate t is always harmonic under these assumptions. 

The coordinate changes to another system of adapted coordinates {t', ip'} pre¬ 
serving the regularity in the axis of symmetry and with closed compact orbits of 
periodicity 27r for the axial Killing vector are t = at' and ip = p>' + bt', where 
a and b are constants. If at spacelike infinity the metric tends to the Minkowski 
metric then a = 1 and this implies for the infinitesimal change of Section 0] that 
= 0,C°(a^) = 0, and taking also into account the independence of the metric 
on time coordinate we have ^{x^) and C''{x’^), i.e. do not depend on time coordinate. 
Furthermore the change ip = (p' + bt' do not maintain the harmonicity condition, i.e. 
if the Cartesian coordinates x, y, z associated to the spherical r, 9, p are harmonic 
the Cartesian ones associated to r, 6 , p' are not harmonic. 

Consequently, splitting time and space components in the equations fl32)) and 
fl3T)l we have 

First order 


0 

0 

0 

0 

( 38 ) 

hojix'^) 

= hoj{x^) 

( 39 ) 

hij{x^) 

= hij{x^) - diij{x^) - dj^i{x^) 

( 40 ) 


the metric components /iqo y hoj of the metric deviation are invariants to first order. 


Second order 


Qooix'") = qoo{x^) - i\x^)dihoQ{x^) (41) 

qojix'^) = qojix’^) - i\x^)dihoj{x^) - hoidji\x^) (42) 

gjj(x^) = qij{x^) - diCj{x^) - djCiix’^) - £{$,)hij{x^) - diC\x’")djCi{x’^) (43) 

With respect to the “harmonics changes” fl36p . the stationarity and axial sym¬ 
metry with a Papapetrou lead to the following equations 


[> 

o 

(44) 

AC(x) = h^Hx)djke{x) 

(45) 

i.e. the functions must be three independent harmonic functions. On the other 
hand, the functions C* must be the sum of a harmonic function and a particular 
solution of equation fl45D. 

Finally, the harmonic condition and Einstein’s equations of Section [3] with the 
splitting of time and space components become: 

First order 

— 27oo * 

(46) 

AA™ = -2T,f ; = 0 

(47) 

a;4‘> = -2Tf >; 

(48) 

where h := S^^hij. 

Order n > 2 


A_ _27"0) I 

AAUqo — -I- ZJVqq 

(49) 

Ahg) = -2T^f + 2N^f - 

(50) 

AhS;) = -2%f + 2 - 

ISiC (51) 


6 Vacuum exterior solution in harmonic 
coordinates 

Let’s assume that we have a compact stationary gravitational source with axial 
symmetry and we would like to study the exterior 2-postminkowskian solution [1]. 
First of all we need to solve equations fl46D . fl47)) and fl48|) with = 0, and then 
to solve the equations (53), (54) and (55) for n = 2 and T/lu = 0. Henceforth we 
will use the notations of (28) and (29). 
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6.1 First order in A 

• For the scalar term hoo the solution of fl46D is the solution of the Laplace equation 
in the spherical coordinates, introduced in Section [5l The solution is 

°° A/f 

hoo(,r,9) = ‘^'^-^PnicosO) (52) 

where Mn are arbitrary constants, representing the static multipole moments of 
Thorne PE], and Pn{x) is the Legendre polynomial of order n. 

• For the vectorial hoj terms the solution of flTTl) is 

OO j 

hoj{r, 0) = 2 ^ ^P,,^(cos e)e^j (53) 

72=1 

where are arbitrary constants, representing the dynamic multipole moments of 
Thorne, := (— sin (p, cos (p, 0) is the tangent vector associated to the (p coordinate 
and Pnix) is the Legendre function of the hrst kind. 

• For the tensorial hij terms the solution of fHS]) is 

hij = hooSij + diWj + djWi (54) 

where 

°° O °° D 

•= ^^nicos9)e,j + ^P„'(cos 9)epj (55) 

72=0 72=1 

and Qn and Rn are arbitrary constants, Cp := (cos </?, sin 9 ?, 0 ) is the tangent vector 
associated to the p cylindrical coordinate and := ( 0 , 0 , 1 ) is the tangent vector 
associated to the 2 ; cylindrical coordinate. 

It is clear that the solution of the homogeneous part of the tensor equation fl5TD 
is pure gauge [see (14(11) ] . That is, it could be interpreted as a first order inhnitesimal 
change of coordinates with ^*(^) furthermore it would be a “harmonic 

change” because it is easily verified that Ate* = 0. On the other hand the complete 
equation, that contain the derivative of the scalar equation solution, has the simple 
particular solution hoo(5ip 

Now assuming that the hrst order solution contains this gauge dependence in 
the tensor term we would like to know how this gauge terms, Wi, propagate to the 
second order term of the deviation, Qag. To simplify the resulting equations we are 
going to write down only the terms containing the hrst order gauge. 

In particular, we would like to verify that the resulting terms coincide with those 
that come from a inhnitesimal change of coordinates up to the second order. 

6.2 Second order in A 

The second order equations to solve are (I4^ . (l50|) and (15TD with = 0 and n = 2 , 
then we need to use equation (1231) where only the deviation to hrst order in A 
appears. 
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Let us assume the gauge dependence in the hrst order deviation given by equa¬ 
tions fl5^ . fl5^ and fIMl) . i.e. only the terms depend on the hrst order gauge. 
From fl26l) we have then in only the terms k^^ depend on the 

hrst order gauge. With this gauge dependence let us compute the terms F^) 

^ 0,00 ~ 0 ) ^o,oj ~ 2^j^oo^ = — {dihoj + djhoi^ ( 56 ) 

^fc,oo “ ~2^khoo: = -{djhok — dkhoj) (57) 

~ 'T^i^kjdihQQ -\- SkidjhQQ — Sijdkhooj + dijWk (58) 

Only the terms Fjand depend on the gauge, now we can write the gauge 
dependent terms splitting the space and time components of and in fl23l) 
and with the sign ~ we mean that we omit the terms with no gauge dependence. 

F^(2) ^ ^ + d^w’^) {dkhoi + dihok) (59) 

pd) ^ ~ -{diWk + dkWi)d^hoo + diW^dkhoo - 2d^w^dijWk (60) 

where we have taken into account that Arc® = 0 

- riprF ~ d‘wia„h„, (ei) 

I^Qj — K Oxi- p Qj J- O/i jX — 

{d’^w^ + d'-w^) {dkjhoi - dkihoj) - ^ [d^h^ - d’^h’-QjdjkWi (62) 


Njf = k^^dxT^;}^ - fS;)^fW" - 1 {k^^d,,hxp + d,hx,d,k^^) ^ 

Sijdkihood'^w^ + dijhoodkW^ - d(^ikhoo{dj)W^ + d^Wj)) 

-h(9^/ioo {dkiiWj) - dijWk) d(ihoodj)kW^ + dikWidjW^ (63) 

where as usual {i j) mean simmetryzation on the index i,j. 

Once we have determined the gauge dependent terms of the right hand side of 
the equations we are going to hnd the second order solution. As we did before let 
us begin solving the scalar equation 

Agoo = ~ 2d'^wWijhoo (64) 

The solution would be equal to the homogeneous solution, formally identical to 
the hrst order more a particular solution of the inhomogeneous equation. But, the 
second exact right hand side of the equation is the sum of the part that does not 
contain the vectors plus the part that contains them. Then it turns out then that 
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the particular solution can be written as the sum of a canonical particular solution 
plus a pure gauge particular solution. Therefore, for our purposes, we need to find 
out a particular solution of the differential equation: 

Aqoo = 2d'"wWijhoo (65) 

taking into account that Hqq and are solution of the Laplace equation then 

A{w"dihoo) = 2&>w^djihoQ ( 66 ) 

Thus, the particular solution w^dihoo has exactly the structure required by the prop¬ 
agation of the hrst order gauge to second order fITT]) . remember that iVj := 

Let us now solve the vector equations 

Aqoj = 2N^f - ~ 2d’^w^dkihoj + 2djkW^d’^hoi 

d^qoj = ^ d^w^{dkhoi + dihok) (67) 

like in the scalar equation the homogeneous solution is formally the same as the first 
order one and taking into account that Awk = 0 and d^hoj = 0 we obtain that a 
pure gauge particular solution is 

w’^dkhoj + hokdjW^ ( 68 ) 

i.e., the solution required by the propagation of the hrst order gauge to second order 
as is given in equation (IT2D . 

The tensor equations are 

Aqij = 2N^f - diPj'^'’- djPp\ Qij - ^diq =-^diqoo - P-^^ (69) 

where 

2Nlf - diPf^ - djP^ ~ 2d'"w^dkihoo6ij + 2d’"hoo{dkiiVj + dkjWi) 

+ 2dikWidjw'' + 2[d^^WidjkWi + d^^WjdikWi) + 2d^w’-{dkiiWj + dkijWi) ( 70 ) 
and using the result of the scalar equation we have 

-^^i<7oo - - - ^diw'^dkhoo - ^w'^dkihoQ + 2d^w'‘dkiWi 

+ d'"hoo{dkWi + diWk) - dkW^dihoo (71) 

Following the same reasoning we have done above, the particular solution depending 
on the hrst order gauge term of the diherential system should be the gauge term 
transferred from the hrst to the second order flT2D . i.e. 

qij ~ w^dkhij + hkjdiw’" hikdjw’^ - diw’^djWk (72) 

where 

hij ho^dij “1“ 3iWj -\- djWi (78) 

and it is easy to see that is a particular solution of equation (16911 . let us be think 
that this is only the particular solution depending on the hrst order gauge, the terms 
depending on the rest of the hrst order solution that appear in P and N and the 
terms depending on the gauge of second order are not computed. 
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6.3 The possible validity of the result to the third order 

If we replace fl73|) in fl72|) . after a short calculation we obtain 

q*j = w’^dkhoo 6ij + hoo{diWj + djWi) + diWkdjw’^ 

di{w''dkWj) + dj{w'^dkWi) (74) 

Therefore the exact solution of fl69ll can be written as follows 

Qij = Qij + w’^dkhoo 6ij + hoo{diWj + djWi) + diWkdjw’^ 

+ di {wf' + w^dkWj^ + dj ^ + w'^dkW^ (75) 

where g))- represent a particular solution of (1^ with tc* = 0 and where the sum 
diwf^ + djwl‘^'^ is the general solution of the homogeneous system, therefore Ate® = 
0 , i.e. with formally identical to (jS^D- 

The second line of fl7T|) could be interpreted as the main part of an inhnitesimal 
change of coordinates to second order. It is also easy to see that Q = + w^dkii 

is the general solution of the second equation of fl45|) with = —ta* and Q = w\ , 
that is it would be an harmonic coordinates change. This suggests that the results 
concerning to the gauge transferred would also be valid to the third order. 

7 Interior solution in harmonic coordinates for 
a rigidly rotating perfect fluid with a linear 
equation of state 

We are going to use the solution obtained in j^]. 

From the metric equation fl371) and 

Cri dx^ = dr, CQi dx^ = rdO, e^i dx* = rsindd(p. (76) 

we have 

ds^ = -fttdf+‘^ltipe^jdtdx^ + {^rrerierj+‘2^1reerieej+^eeeeieej+lipipe^ie^j)dx"dx^ (77) 
that is 

9ij ■ '^rr^ri^rj T T '^OO^di^dj T i 90j ■ 

We assume that the source of the gravitational held is a perfect huid, 

Ta0 ■.= To,0- = (/i + p)UaU0 + - p)9a0 (78) 

whose density p and pressure p depend only on r and 9 coordinates. We also assume 
that the huid has no convective motion, so its velocity lies on the plane spanned by 
the two Killing vectors 

u = ii{dt + ujd^) = + ursine e^^di) (79) 
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from now on we are going to assume that the fluid rotates rigidly co = constant and 

= [- (yti + 2u 'Jtip r sin 6* + sin^ 9) ] ^ (80) 

is a normalization factor. In order to obtain a postminkowskian development for 
the metric we need to have the density fi of order A i.e. fi = \fl. 

In the postminkowskian approximation we have 

900 = —1 + Ah-oo + A^goO) 9oj = <^(Ahoj + A^goj) > 9ij = ^ij + + A^g^j (81) 

this implies that 

'ytt f ^ Qtt'i 'ytif “^(Ah^i^ A 9 t<p^^ ^ A Q^p^p 

(82) 

7rr = 1 + Xh^r + A^g^r, 7r0 = Xh^-e + X^Qrd, 'Jeg = 1 + Xhgg + X'^ qgg 


Let us consider Euler equations for the fluid (or, what is equivalent, the energy- 
momentum tensor conservation law) [TT] 


daP = (h + P)da In V' {a,b, - ■ ■ = r ,9). 
we are going to use the following linear equation of state (EOS) 

fi + {l-n)p = po 

with this EOS the Euler equations are integrated easily, to give 




n 




P=-U^ -1 






n 


P=ziiri-1) ^ +1 






(83) 


(84) 


(85) 


where is the value of the normalization factor ijj on the surface of zero pressure, 
which in turn leads to the following implicit equation for the matching surface 


p = 0 




( 86 ) 


with this dependence on the parameter, the lowest order in A of this equation give 
to us cylindrical surfaces instead of spherical deformed ones unless the rotation 
parameter depends on A, we need oc A at least then to the lowest order in A 

■0 ~ 1 + -1-(h^r^sin^ 0 ) (87) 

where we have defined := Ao)^ and for the constant 

V's := 1 + A5 (88) 

= Xu"^ {u oc n, slow rotation parameter) so that this approach is consistent with 
the post spherical approach carried out in previous articles (see details in references 
d] and [2]). This implies that some quantities, including gQj, have a series expansion 
in half powers of A. However, this situation does not affect the present work, in 
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which we only deal with the postminkowskian approach, so we are just to consider 
the integer part of the orders of the series in the computation. 

Then the approximate pressure is 


p — SS* + htt + sin^ 6*) 


(89) 


and for the velocity field we have 


Mo ~ — 1 + sin^ 0) 


Uj ^ <; r sin 6^ + A 


htip + -r sin 6>(htt + 2h<^<^ + ui^r^ sin^ 9) 


(90) 


■'w 


now we can compute the tensor energy momentum to second order in A and we 
obtain 

Too ~ + ^A^/i( - 2(n + 2)5” + nhoo + {n + 6)a;^r^ sin^ 6 ) 


Toj ~ sin Oe^pj — 

hofcC^ + r sin 6'( - 2nS + nhoo + 2hfcte^e(, + (n + 2)a;^r^ sin^h) 




(91) 




+ tA^p (2 — ri) (25” — h-oo — sin^h)(jy + 2/iy + sin^ de^ie^j 


7.1 First order in A 

The Einstein’s equations to first order are (1^ . fl47l) and fl48l) . the scalar equation 

Ahoo = -2%^o^ = -p (92) 


has as solution, well behaved in r = 0 given by 

^ OO 

hoo{r, 6) = --jl + 2 ^ Pn(cos 9) 

n=0 

The vector equations are 

Ahofe = —2TqI'^ = +2fb ursm9 e<^fc 

d‘‘hat = 0 


(93) 


(94) 


and their solution is 


hok{r, 9) 


:fiur sm 


h + 2^T„r”P^(cosh) 

n=l 


^ipk 


(95) 
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as for the exterior solution there are not gauge dependent terms in the hrst order 
solution for hoo and hoj. 

The tensor equations are 



Ahij = = -ftSij 

d'^hkj -^djh = -^djhoo 

(96) 

and their solution is 

hij = hoo 5ij + diWj + djWi 

(97) 

where now 

OO 

OO 

^r^Pn{cos9) Czj + 'Y Rnr'^Pnis^osO) Cpj 

(98) 


n=l n=l 


As in the exterior problem the solution of the homogeneous part of this tensor equa¬ 
tion is pure gauge, see (HUj) with := —Wj . And the complete equation that contain 
the derivative of the scalar equation solution has the simple particular solution hoo^ij 
too; then we can use the same expressions flSU]) . (IMD . (1^ and (1^ we have 

used above for the exterior solution. 

Now also we want to determine how the terms “pure gauge” of the hrst order 
solution are transmitted to the second order solution and in particular verify that 
the resulting terms coincide with those coming from an inhnitesimal second order 
change of coordinates. 

7.2 Second order in A 

We can now go to second order i.e. to solve fl4^ . flHHH and flHTD with T given by 
equations d^T]) when n = 2. Do not confuse this n related with the order of the 
approximation we perform with the n appearing in the EOS. As in the exterior 
problem we are going only to look for the terms in the solution corresponding to the 
propagation of the hrst order gauge. 

The scalar equation is 

Agoo = 2ivS^ - 2rolf^ ~ 2d^w^d^jhoo (99) 

/o'! 

the term Tqq do not contain the gauge vectors tc* the terms to solve have the same 
aspect that in the exterior case. And we have 

A(tD*9ihoo) = Ate* dihoQ + 2d^w'^ dj{dihQQ) + tv^A{dihQQ) = 28^dphoQ (100) 

where we have taking into account that 

Aw* = 0 , A(aihoo) = 9i(Ahoo) = 0 (101) 

Then, the particular solution for the terms depending on the hrst order gauge has 
the aspect required for a gauge change to second order fl34)l . 
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For the vector equations the term only is gauge depending in the hij term, 
that is 

Tof ~ -2ilujr sin e^j (102) 

the equations to find a particular solution of the gauge depending part are 

A<;oj = 2]vg> - ajPf' - 2r„f 

~ 2d^w^dkihoj + 2djkW^d^hoi + Ajlur sin 9 {dkWie^^e^^)e^j (103) 
d^qoj = ^ d^w^{dkhoi + dihok) 

We know from the second order gauge infinitesimal transformation fl34j) that the 
expression w^dkhoj + hokdjW^ should be the right answer, let’s try this particular 
solution 

A{w^dkhoj + hokdjW^) =2d^w’^ dikhoj + 2d^hok dijw^ 


+w^dkAhoj + Ahokdjw’^ (104) 

and using the first order vector equation fIM)) and the expression for Uk we get 

CX) 

w^dkAh^j + AhokdjW^ = 4:fluj (cos 6*) e^pj 

n=l 

and the last term in the expression for Aq^j using 

OO 

rsine{dkWie^^e'‘^) = ^ P„r”P„^(cos 6^) (105) 

n=l 

Now if we put this particular solution in the harmonicity condition 

cF' [w^dkhoj + hokdjW^) = {djhok + Skh^j) (106) 


then the gauge propagation gives the particular solution to second order for the 
vector equations too. 

To finish with let us solve the following tensor equations 

A,., = 2/vf - 9,pf» - a,P^’ - 27-ft (Pqti - - ff > (107) 

the only difference with the equations we have solved for the exterior problem is the 
term 

(108) 

and the remaining terms are determined from fl63|) and fl60l) and the solution of the 
scalar equation (1^ . 

Then as we did above let’s try with the particular solution 

q*j = w^dkhij + hkjdiW^ + hikdjW^ - diW^ djWk (109) 

where 

hij = hoo6ij + diWj + djWi ( 110 ) 

As in the exterior vacuum problem a direct calculation shows that this is indeed 
solution. 
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8 Conclusions 


When we solve the Einstein equations in harmonic coordinates to each order we 
obtain solutions that depend on the development parameter and some constants, 
we can distinguish to types of constants multipole moments and the so called gauge 
constants, giving all solutions but keeping them within the family of harmonic coor¬ 
dinates. These gauge constants are essential if we match the interior solution to the 
exterior solution to get a global solution but with this procedure we obtain messy 
solutions when we compute the nonlinear terms in the development constant A. 

It seems natural to think that an infinitesimal change of harmonic coordinates 
depending on the same constant A allows us to transform the metric and simplify 
the solutions so that we can manipulate solutions that have not gauge constants and 
obtain the general solution by performing an inhnitesimal coordinate transformation. 

This paper demonstrates that this happens to at least second order in the con¬ 
stant A. 

It is to say, we can find a hrst order solution with all the gauge constants then 
with this solution we can generate the second order solution which will contain 
quadratic terms in the hrst order gauge constants and on the multipolar hrst order 
constants and linear terms on the second order gauge constants, in n. 0 we do this 
following the procedure sketched in the appendix, let’s call this solution g or we can 
hnd the hrst order solution with the hrst order gauge constants equal to zero and 
then with this solution generate the second order solution that only contain multipole 
terms. Then we can transform this solution without gauge terms by a second order 
inhnitesimal coordinate transformation, with the condition that transforms from 
one set of harmonic coordinates to another set of harmonic coordinates, let’s 
call this solution g then with a redehnition of the gauge constants and of the second 
order multipole moments we can made to coincide the solutions g and ^0 

This last procedure, without gauge constants, to hnd the solution let us to get 
simpler solutions and saves computing time, then if we we need to match them with 
another interior or exterior solution we must perform a second order inhnitesimal 
harmonic coordinates transformation on the solution without gauge constants to 
obtain a solution with all the constants needed to do the Lichnerowicz matching of 
solutions. 
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